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I. WHAT IS FREQUENCY RESPONSE? 

1. Recall Sinusoidal Steady State Analysis 
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EXAMPLE 2. A bandpass (BP) transfer function is one that 
passes frequencies in a (narrow) frequency band and 
eliminates/attenuates-significantly those freqs outside the band. 



K cos(cot + 0 ) 

Input ► 

Sinusoid 



HO Oft) 



M(Cfl) COS(G0t + <P(G>)) 
► 



Steady State 
Output 
Sinusoid 



(a) M{(0) = Kx \H(jco)\ - to -dependent gain 

(b) (p(co) = ZH(jco) + 6 - co -dependent phase 

2. Definition: the frequency response of a circuit/system 
is H(jco) = H(s)] s _ i.e., H(s) evaluated along the imaginary axis. 

NOTE: for each co, H( jco) is a complex number: 
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(a) ZH(jco) isthephase response, and 

(b) \H(jco)\ is the magnitude response /GAIN. 

5 + 10 



EXAMPLE 1. Plot the Phase and Magnitude of H(s) = 



5 + 400 



S = JCO 



\H(jco)\ 



ZH(jco) 



j0 



0.025 



no 



10V2 



/10-+400 Z 



10V2 
'400 



, 10+./10 

"400 + /10 
= 45° -0 = 45" 



7400 



400 



400V2 V2 



, 10 + 7400 

"400 + j'400 

= 90" - 45" = 45" 




Frequency r/s 

n = [1 10]; % numerator coefficients 

d = [1 400]; % denom coefficients 

w=0. 1:0. 1:2000; % Frequency range 

h = freqs (n,d,w) ; % Freq response calculation 

figure (1) ;plot(w,abs(h) ) ; grid ;xlabel ( 'Frequency r/s') ;ylabel ( 'Magnitude response') 
figure(2) ;plot (w, angle (h) *180/pi) ;grid ;xlabel(' Frequency r/s') ; ylabel(' Phase response') 



(i) 2 nd order BP transfer function (blue curve): 
0.255 0.255 



5 2 +0.255 + l (5 + 0.125) 2 +(0.992) 2 



(ii) 4 th order BP transfer function (green curve): 

H (s) = 0-Q625. 2 

2 ^ s 4 + 0.35355s 3 + 2.0625s 2 + 0.35355s + 1 



V0.0625s 



Vo .0625 s 



(s + 0.0962 ) 2 + (1 .0884) 2 (s + 0.08058) 2 + (0.91 164) 2 



Remark: s - jco - 0 and s - jco = °° are the two most important 
frequencies: 

(a) //j (0) = H { (joo) = 0 (magnitude is zero at dc and infinity) 

(b) ff 2 (0) = # 2 (» = 0 

(iii) Plot \H(jco)\ using the following MATLAB code: 

»nl = 0.25*[10]; % numerator coefficients HI 

»dl = [1 0.25 1]; % denom coefficients H2 

»n2 = 0.0625*[1 0 0]; % numerator coefficients H2 

»d2= LI 3.5355e-01 2.0625 3.5355e-01 1J; 

»w=0 .2:0.005:2; % Frequency range 

»hl = freqs(nl,dl,w); % Freq response calculation 

»h2 = freqs(n2,d2,w); 

»plot(w,abs(hl),w,abs(h2)) 

»grid 

»xlabel('Frequency r/s') 
»ylabel('Magnitude response') 
»gtext('2nd Order BP') 
»gtext('4th Order BP) 
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(b) Qualitative Analysis Using Pole-zero plot. Imagine, as you 
walk along the imaginary axis, what happens to your height 
as you near a pole, a double pole, and then as you near a 
zero. 




(c) The Idea of Frequency Scaling: 

»Kf = 1000; % Frequency Scale Factor 

»nl =0.25*[1/Kf0]; % s s/Kf in num coefficients of HI 

»dl=[l/Kf A 2 0.25/Kf 11; % s^ s/Kf in denom coefficients of HI 

»n2 = 0.0625*[ l/Kf A 2 0 0]; % s s/Kf in num coefficients of H2 

»d2 = [1/KfM 3.5355e-01/Kf A 3 2.0625/Kf A 2 3.5355e-01/Kf I]; 

»w=0.2*Kf: 1 :2*Kf; % Freq range is scaled up proportionately 

»hl = freqs(nl ,dl,w); 

»h2 = freqs(n2,d2,w); 

»plot(w,abs(hl ),w,abs(h2)) 

»grid 

»xlabel('Frequency r/s') 
»ylabel('Magnitude response') 
»gtext('2nd Order BP') 
»gtext('4th Order BP') 
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EXAMPLE 3. Magnitude response \H( jco)\ of three (so called) 
normalized low pass (LP) Butterworth filter transfer functions: 

(a) First Order Normalized LP Butterworth 




H(s)-- 



s + \ 



(b) 2 nd Order Normalized Butterworth (how many 00 zeros) 
1 



H(s)=- 



LC 



1 



2 + ^ + _L s 2 + j2s + l 
L ' LC 




ont 



(c) 3 rd Order Normalized Butterworth (how many 00 zeros) 
1 



H(s)~- 



r+2r +2s+i 




(d) MATLAB Code for frequency response calculation given H(s): 

»w = logspace(-2,2,800); 

»nl = l;dl = [1 I]; 

»n2= l;d2 = [1 sqrt(2) 1]; 

»n3 = 1; d3 = [1 2 2 1]; 

»hl = freqs(nl ,dl ,w); 

»h2 = freqs(n2,d2,w); 

»h3 = frcqs(n3,d3,w); 

»semilogx( w ,abs(h 1 ) ,w ,abs(h2) ,w ,abs(h3 )) 

»grid 

»xlabel('Normalized rad frequency') 
»ylabcl('Magnitude Response') 



(e) Magnitude Response Plots (Note transitions of gain from 1 to 
near 0) 




Normalized rad frequency 
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(f) Pole-Zero Plot of 2 nd and 3 rd Order (Normalized) Butterworth 
Filters— Notes: How many infinite zeros? Poles lie on the unit 
circle in left half complex plane. Are such filters stable?? 
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EXAMPLE 5. Frequency response of the (High Pass-HP) RC 
circuit; Capacitor has hight impedance for low frequencies and low 
impedance for high frequencies. 




3rd Order 
Poles 



2nd Order 
Poles 



EXAMPLE 4. Plot of \H(jo))\ of three low pass Butterworth 
transfer functions of frequency scaled by K j ■ = 1000 : 

I 



H new( s )~ H 



S 



old 



v K f , 

\ J J 



— v> 



(a) Transfer function: H(s) = - 



5 + 100 



by V-division. 



C 



(i) One pole at s = -100 ; one zero at s = 0 . 

-— 1(- 



in 



6 



0.01F 



ont 



(b) S = j<D,H(j<D) = 



jcoCR 



. (O 

J Too 



* + _L jcoCR+x jjL+i 

jG)C 100 



Step 2. Asymptotic Behavior: Important frequencies are: 
0,oo, and ???? 



H new (s) = H 



( \ 
s 



old 



H neA s )~ H 



V K f , 

V J J 



f \ 

s 



old 



f— V' 

Uooo J uooo ) 



f— T**f— IV— 1 

uooo J uooo ) uooo J 



+ 1 



Note: Same Plots as before with cut-off frequency at 1000. 




10 3 

Normalized rad frequency 



6l)->100 



\H(j(0)\ 



R + - 



jcoC 



-> 1 



\HU®)\ = 



jcoCR 



jct)CR+\ 



->0 



\H(jco)\ = 







J \o6 
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/ + 1 




' 100 





ZH(jco)^0 



ZH(jo))^ 90° 



ZH(j<0)-> 45° 



Step 3. Plots from Matlab: 



»w = logspace(0,4,500); 

»H = freqs([l 0]J1 100],w); 

»semilogx(w,abs(H)) 

»grid 

»xlabel( 'Frequency in rad/sec') 
»ylabel( 'Magnitude') 



High Peas Response 



1(f 

Frequency in radf sec 
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»semilogx(w,180*angle(H)/pi) 
»grid 

»xlabel('Frequency rad/sec') 
»ylabel('Angle in degrees') 
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Step 2: Important Frequencies and Asymptotic Behavior 




EXAMPLE 6. Frequency response of the (Band Pass-BP) RLC 
circuit 




L-0.1H c=10" 3 F 



1 o o> y om 



Step 1: Calculation of transfer function: By V-division, 

R 



H(s) = 



R 



VinW 



or 1 2 R 1 

R+Ls + — s +—S+ 

Cs L LC 



NOTE: a finite zero at s = 0 and one infinite zero; 
two finite poles at s = -50 ± y'86.6 



(0 

(rad/sec) 


.0) 

H{j(0) = -: £ = 4 J , 

LC L 




\H(jQ))\ZH(jco)^0Z90 0 




\H(jco)\ZH(jco)^0Z0° 


0) -> 100 


\H(j(O)\ZHUco)^\Z0° 


77 


\H(jG))\ZH(j(D)^-^Z45° 


77 


\H(jco)\ZH(jo))^-^Z- 45° 



Hint: for what values of co are the magnitudes of the real 
and imaginary parts of the denominator equal? 

Step 3 . Plots from Matlab: 

w = logspace((),4,l(X)0); 

R= 10;L = 0.1;C = le-3; 

n = [R/L0]; 

d = [l R/L 1/(L*C)]; 

h ■ freqs(n,d,w); 

semilogx(\v%abs(h)) 

grid 

x label ( 'Frequency in rad/s') 
y label ( 'Magnitude Respon se' ) 
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NOTE: a finite zero at s = 0 and one infinite zero; 
two finite poles at s = -50 ±786.6 
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»semilogx( w Mingle(h)* 1 80/pi) 
»grid 

»xlabel( 'Frequency in rad/s') 
»ylabel( 'Phase Response' ) 




Frequency in rad/s 

EXAMPLE 7. Frequency response of the (Band Reject-BR) RLC 

circuit: Passes most frequencies with little or no attenuation and 
rejects a small band of frequencies. 

L-0.1 H 



+ \ C= 10" J F < Y, 



LQ > 



out 



1 O) 
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10 
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StepI: Calculation of tramsfer Function: By V-division, 



H(s) = 



2 1 
S + — 
LC 



1 2 1 1 

/? + r- .r+ S+- 



J_ RC LC 

Ls 



±7 



Remark: Zeros at ^— — = ±y*100 , frequencies around ± j 100 are 



mostly gone: H(jco) = 



-6)~ 



1 2 J® 

co +- — 

LC RC 



w = logspace (-1, 5, 1000) ; 
R = 1; L = 0.1; C = le-3; 
n = [1 0 1/ (L*C) ] ; 
d = [1 1/ (R*C) 1/ (L*C) ] ; 
h = f reqs (n, d, w) ; 
figure (1) 

semilogx (w, abs (h) ) 
grid 

xlabel (' Frequency in rad/s') 
ylabel ( 'Magnitude Response') 
figure (2 ) 

semilogx (w, angle (h) *180/pi) 
grid 

xlabel (' Frequency in rad/s') 
ylabel (' Phase Response 1 ) 



